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ABSTRACT 

Topic of this article are tomographic measurements of the integrated Sachs-Wolfe effect 
with specifically designed, orthogonal polynomials which project out statistically independent 
modes of the galaxy distribution. The polynomials are contructed using the Gram-Schmidt or- 
thogonalisation method. To quantify the power of the iSW-effect in contraining cosmological 
parameters we perfom a combined Fisher matrix analysis for the iSW-, galaxy- and cross- 
spectra for wCDM cosmologies using the survey characteristics of PLANCK and EUCLID. 
The signal to noise ratio has also been studied for other contemporary galaxy surveys, such 
as SDSS, NVSS and 2MASS. For the cross-spectra our tomographic method provides a 16% 
increase in the signal to noise ratio and an improvement of up to 30% in conditional errors 
on parameters. Including all spectra, the marginalised errors approach an inverse square-root 
dependence with increasing cumulative polynomial order which underlines the statistical in- 
dependence of the weighted signal spectra. 

Key words: cosmology: large-scale structure, integrated Sachs-Wolfe effect, methods: ana- 
lytical 



1 INTRODUCTION 

The integrated Sachs-Wolfe (iSW) effect is one of the secondary 
anisotropics of the cosmic microwave background (CMB). Time- 
evolving gravitational potentials in the lar ge-scale structure gen - 
erate temperature fluctuations in the CMB jSachs & Wolfdll967h . 
The iSW-effect is a valuable tool for investigating dark energy and 
non-standard cosmolo gies since it is sensitive to fluids with non- 
zero equation of state dCrittenden & Turoklll996t) . For this reason 
its detectio n is of particular relevance fo r cosmology and the nature 
of gravity dLue et al.ll20ol : IZhandl2006l) even if its signal strength 
is very low. 

Since the iSW-effect is generated in time-evolving poten- 
tial wells for photons on their way from the last scattering sur- 
face to us, it will be strongly correlated with the galaxy den- 
sity field. Therefore, the cross-spectrum will provide valuable 
additional cosmological information. The i SW effect has been 
measured in such cross-correlation studies dBoughn et al. 1998; 



Boughn & Crittenden! |2004|; IVielva et alj |2006| ; iMcEwen et al.l 



20071 : Giannantonio et al. 2008T). However, due to the line of sight 



integration, a detailed distance resolution of the processes can not 
be withdrawn from these spectra. 

A former approach correlated large scale structure observa- 
tions from various survey with the CMB anisotropies to study the 
iSW-effect as a function of redshift and to for mulate a reliabl e like- 
lihood f ormulation for paramete r constraints dHo et al.l2008h . Also 
recentlv iFrommert et al.l J2008) presented an optimal method to re- 



* gero .j uergens @ stud. uni-heidelberg . de 



duce the local variance effect and gained 7 per cent in the signal to 
noise ratio for the cross-spectra. 

In this work we aim to formulate a tomographic approach 
with help of an orthogonal set of weighting polynomials, which 
is similar to a former a pplication to weak lensing spectra 
( Scha efer & Heisenberd |201 l|). The orthogonality of the polyno- 
mials will generically lead to a diagonal signal covariance matrix 
and will therefore provide cumulative statistical independent mea- 
surements with increasing polynomial order. 

The article has the following structure: In Section [2] we 
provide introductory information about dark energy cosmologies, 
CDM power spectra and linear structure growth within these cos- 
mologies (Sections 12, l|2,3t . We also introduce a galaxy distribu- 
tion function (Section I2.4t and give a short introduction to the 
iSW-effect (Section 12.5b . The orthogonal polynomials are moti- 
vated and constructed in Section [3~Tl and [3~2l also their most impor- 
tant properties are discussed (Section[33J. In Section|4]we discuss 
how tomography with orthogonal polynomials can improve statis- 
tical constraints on cosmological parameters. After calculating the 
noise contributions (Section r4.lt we perform a Fisher matrix anal- 
ysis (Section [4.2t and discuss signal to noise ratios and statistical 
errors (Section r 4.3|4.4t . The results are summarised in Section[5] 

The reference cosmological model used is a spatially flat 
vvCDM cosmology with Gaussian adiabatic initial perturbations in 
the cold dark matter density field. The specific parameter choices 
are n ra = 0.25, n s = 1, <r 8 = 0.8, D b = 0.04 and Ho = 
100 h km/s/Mpc, with h = 0.72. The dark energy equation of state 
is set tow = -0.9 and the sound speed is equal to the speed of light, 
c s = c. 
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2 COSMOLOGY AND ISW-EFFECT 
2.1 Dark energy cosmologies 

In spatially flat dark energy cosmologies with the matter density 
parameter Q. m , the Hubble function H(a) = dlna/dt is given by 



^ n „-3 , ia n \ -3(l+w) 

—p r = il m a +(l-£2 m )a 



(1) 



with a constant dark energy equation of state parameter w. The 
value w = -1 corresponds to the cosmological constant A. The 
relation between comoving distance X and scale factor a is given 
by 



X = c 



I da — , 

J a a 1 H{a) 



(2) 



in units of the Hubble distance^ = c/Hq. 



In contrary to the pressureless dark matter component the 
baryons inside a dark matter halo can loose energy via radiative 
cooling and form stars. Because of this different behaviour, strictly 
speaking, one can not deduce the fractional perturbation An/(n) in 
the mean number density of galaxies (n) from the dark matter over- 
density 6 = Ap/p. In a very simple way, however, the linear relation 
between the two entities, 



An 



Ap 

(P) 



(8) 



is a good appro x imatio n in most cases and was proposed 
by iBardeen et alj dl986|). The bias param eter b can gener- 
ally depend on scale I Lumsden et al.l ll989T) . time dFrvl 1 19961: 
iTegmark & Peebles! 1 19981) as well as the galaxies luminosity and 
morphology. For simplicity we set the galaxy bias to unity through- 
out this paper, b = 1, An established parametrisation of the redshift 
distribution n(z) dz of galaxies is 



2.2 CDM power spectrum 

The linear CDM density power spectrum P(k) describes the fluc- 
tuation amplitude of the Gaussian homogeneous density field S, 
{S(k)S*(k')) = (2n) 3 S D (k - k')P(k), and is given by the ansatz 

P(k) oc k ni T 2 (k), (3) 

with the transfer function T(k). In l ow-fl m cosmologies T (k) is ap- 
proximated with the fit proposed bv lBardeen et al.l dl986l) . 



T(q) 



ln(l +2.34q) 
2.34 q 



X [l + 3.89,7 + (16.1 a) 2 + (5.46a) 3 + (6.71 qf] " 4 (4) 

wher e the wave numb er k = q T is rescaled with the shape parame- 
ter T ([Sugiyama 1995) which assumes corrections due to the baryon 
density Sl b , 



T = Q m /i exp 



-a 1 + 



V2/7 



(5) 



The spectrum P(k) is normalised to the variance tr 8 on the scale 
R = 8 Mpc//t, 



dkk 2 P(k)W 2 (kR) 



(6) 



with a Fourier transformed spherical top hat filter function, W(x) = 
3/i(j)/ x. jt(x) is the spherical Besse l function of the first kind of 
order t dAbramowitz & Stegunlll972h . 



2.3 Structure growth with clustering dark energy 

Linear homogeneous growth of the density field, 6(x,a) = 
D + (a)6(x,a = 1), is described by the g rowth function D + (a) , 
which is the solution to the growth equation l Turner & White! 19971 : 
IWang & Steinhardtll 19981: iLinder & Jenkinll2003h . 

d 2 1 / dkiff\ d 3 

— D + (a) + - 3 + — — D + {a) = — ii m (a) D + (a) . (7) 

da 1 a\ din a / da 2 a 1 



2.4 Galaxy distribution 

Galaxies form when strong peaks in the density field decouple from 
the Hubble expansion due to self-gravity. These s o called protoha- 
los approximatel y undergo an elliptical collapse dMo et alj|2007t 
ISheth et alj|200lh . 



with I = * r f!) 
no P \P) 



(9) 



n(z) dz = n — exp - — dz 
\za) [ \zo) _ 

wich was introduced bv lSmailetai]dl995h and will also be used in 
this work. The parameter zo is related to the median redshift of the 
galaxy sample z me d = 1.406 z o if /J = 3/2. Finally, the T-function 
dAbramowitz & Stegunll972h determines the normalisation param- 
eter n Q . 



2.5 The integrated Sachse- Wolfe (iSW) effect 

Due to its expansion our universe had cooled down sufficiently to 
allow the formation of hydrogen atoms at a redshift of z — 1089 
dSpergel et alj2003l) . Fluctuations in the gravitational potential im- 
posed a shift in the decoupled photons which were emitted in 
the (re)combination process (Sachse-Wolfe effect). This primary 
anisotropy can be observed in the cosmic microwave background 
(CMB) in form of temperature fluctuations AT/Tcmb — 10~ s on 
large scales around its mean temperature 7cmb = 2.726 K dFixsenl 
l2009h . 

Besides this, photons are subjected to several other ef- 
fects on their way to us, which lead to secondary anisotropies 
( Agha nim et alj |2008), of which only the most important ones are 
menti oned here: Gravitational lensing dBartelmann & SchneideJ 
1 20011) . Compton-c ollisions with free cluster el ectrons (Sunyaev- 
Zeldovich effect, IZeldovich & Sunvaevl Il980h and with elec- 
trons in uncollapsed structures (Ostriker-Vishniac effect, 
Ostriker & Vishnia & Il986l) and gravitational coupling to lin- 
ear time-evolving po tential wells (integrated Sachs-Wolfe effect, 
Sachs & Wolfe 1967, which will be subject of this work). 

Assuming a completely transparent space, i.e. vanishing op- 
tical depth due to compton scattering T opt (q) = 0, the temperature 
fluctuations t{Q) generat ed by the iSW-effect can be expressed by 
the line of sight integral dSachs & Woifelll967l) 



T(0): 



AT., 



c Jo 



d X a 2 H(a)—<t>(0x,x) 
oa 



(10) 



reaching out to the limit of Newtonian gravity. Using the Poisson 
equation we can write this integral in terms of the dimensionless 
potential tp = ^lx\ = ^ ' ^lx\ f rom the density field 8: 



r(0) 



3£1 T 



, d D + 

dxa 2 H(a) — -^<f>(8x,x) 
da a 



(11) 



Heuristically, the effect originates from an unbalance between the 
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photon's blue-shift when entering a time varying potential well and 
the red-shift experienced at the exit. 

The effect vanishes identically in matter dominated universes 
Q m = 1, since then D+/a is a constant. Therefore, a non-zero iSW- 
signal will be an indicator of a cosmological fluid with w + 0. 
After the radiation dominated era it will thus be a valuable tool for 
investigating dark energy cosmologies. 

Since the inverse Laplacian which solves for the potential in 
the Poisson equation introduces a k~ 2 term, small scale fluctuations 
will be quadratically damped. For this reason the iSW-effect pro- 
vides a signal on large scales and will be negligible above ( « 100. 

In order to identify the sources of the effect it is sensible to 
investigate the cross correlation of the iSW amplitude with the line 
of sight projected galaxy density y: 



7(0) = b 



f 

Jo 



4f 



(12) 



We obtain the dimensionless observables y and r from a line of 
sight integration of the two dimensionless source fields 5 and <p 
weighted by functions which carry units of inverse length. 

If one is interested in rather small scales one can approximate 
the sphere locally as beeing plane and perform a Fourier transform 



y(t) 



J d 2 7 (0)e-^- fl ) 



Clearly, there is no directional dependence, y(t) 
can define the spectrum C yy (l): 

<7(£)7*(0> = (2x) 2 6 D ({-£')C yy (i) 



(13) 

yif), and one 
(14) 



The observable r can be transformed in analogous way. With the 
two weighting functions 

W y ix) = n(z) — D + (z) 
c 



c aa a 



(15) 



we can now derive the spectra dLimbeifl953h . 

C XH dy , 
Jo X 

CryV) = f X " % W T (x) W y ix) P S4 ,(k = tlx) 
Jo X 



C TT ({) 



Jo X 



(x)PU k = i lx) 



(16) 



The power spectra can be related to the density power spectrum: 



Pss(k) 



Ps S (k) 



(17) 



The multiplication factors k~ 2 and £r 4 tilt the spectra to smaller 
values for increasing mutipole order I and show once again the 
iSW-effect to be a large scale phenomenon. 



3 TOMOGRAPHY WITH ORTHOGONAL 
POLYNOMIALS 

3.1 Motivation 

Measurements of the iSW-effect provides integrated information 
about the structure formation history of our universe since the last 
scattering surface. Cross-correlation with the galaxy density field 
increases the signal to noise ratio significantly and the spectrum 
is noiseless due to uncorrelated noise in the CMB and the density 
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Figure 1. Orthogonal polynomials Pi(x), < = 0...4, as a function of comov- 
ing distance X- The lowest order polynomial is shown in blue, the highest 
order in green. The construction was performed with the Gram-Schmidt 
algorithm at multipole order I = 100. 



field. However, due to the fact that both the cross-correlation spec- 
trum and the galaxy spectrum are line of sight integrated quantities, 
non-linear effects of parameters on the signals could be averaged 
out and valuable tomographical information would be lost. 

Tomographical methods split up the signal from different dis- 
tances and are therefore able to increase the signal to noise ratio and 
the sensitvity with respect to cosmological parameters. In case of 
the galaxy spectra this implies that additional covariances between 
the different spectra have to be taken into account. 

For a direct tomograpphy in the line of sight integration of 
the iSW signal the knowledge of the large scale structure poten- 
tial would be necessary. A reconstruction of the potential from the 
galaxy, however, would not reach the required accuracy. 

To circumvent this issue we perform tomography in the galaxy 
signal and cross-correlate these with the iSW signal. In the course 
of this we are able to withdraw tomographical information also 
from the iSW signal. We use specifically designed polynomials 
for a distance weighting of the galaxy distribution. Defining the 
weighted galaxy covariances as a scalar product of the polynomials 
will lead to statistically independent galaxy spectra once the poly- 
nomials are orthogonalised. This nonlocal binning of the galaxies 
leads to a diagonalisation of the galaxy signal covariance matrix. 
The polynomials can then also be used for tomographical measure- 
ments in the iSW-galaxy cross-correlations. 

3.2 Construction of orthogonal sets of polynomials 

Weighting the given galaxy distribution function nix) = 
n(z)dz/dx = n(z)H(z)lc with a polynomial Piix) modifies the 
galaxy weighting function to 

W?(X) = Piix) W y (x) = Pi(x) "(z) ^ D + (z) . (18) 
For the polynomials p,0t) and pj(x) we require orthogonality 
<Pi,pj> = far (19) 
with respect to the following scalar product for the polynomials: 



J"XH Ay 
^wf(x)w?(x)P(k = eix) 
o X 



(20) 
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polynomial order i+ 1 

Figure 2. Numerical accuracy for the orthogonality relation (pi,Pj) at t = 
20 in logarithmic representation. The accuray imposes a limit on the number 
of included polynomials. 
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Figure 3. Weighted galaxy efficiency function Wj (x)> ' — 0...4, as a func- 
tion of comoving distance at multipole order t = 20. 



The necessary properties for a scalar product are obviously full- 
filled ((pi,pi) > 0, (pi,pi) = «=> p, ■ = and linearity). We use 
the Gram-Schmidt procedure to construct orthogonal polynomials 
out of the family of monomials 



(21) 



where XnoAc sets the position of the node of the first polynomial, 
which is in our case set to the median value of the redshift distri- 
bution. However, a change in XnoAe is completely absorbed in the 
coefficient and has no influence on the polynomials. Starting with 
the zero-order polynomial 



Po(x) = PoOd) = 1 > 

the polynomials are constructed iteratively, 



p,(x) = p',(x)-Yj 



(.P'i-Pj) 



(22) 



(23) 



The procedure has to be performed for every multipole I. The index 
£ of the polynomials Piix) has been omitted for clarity. As one can 
see, the zero-order scalar product is equal to the galaxy spectrum: 



(P0,P0) = Syy(t) 



(24) 



Therefore, the unweighted case is already contained in the first 
weighting function. Finally, we can weight also the tracer density 
modes y{€) themselves with a polynomial Piix) 



f 

Jo 



(25) 



for which a generalized version of the well known expression for 
the covariance holds in case of homogeneous and isotropic random 
fields: 



roughly at the positions where the previous polynomial reaches a 
local maximum or minimum, which intuitively indicates their or- 
thogonality. 

As one can see in Fig.[2]orthogonality is fulfilled until numer- 
ical limitations become significant at a polynomial order of q » 9. 
The inreasing numerical deviations from the orthogonality condi- 
tion ({pt, pj) = for i t j) is due to the iterative method, which cu- 
mulates errors throughout the process. This implies the accuraccy 
to shrink from 1(T 15 for i = to 1CT 3 for i = 8. This is a well known 
disadvantage of the Gram-Schmidt orthogonalisation method, es- 
pecially when dealing with functions as opposed to vectors, since 
there ist larger numerical noise in the evaluation of the scalar prod- 
ucts. However, as we will later see, it is not necessary for our ap- 
plication to go to even higher orders. 

In Fig.[3]the weighted galaxy efficiency functions Wy (x) are 
depicted, which are modified by the polynomials Piix) at an angular 
scale of t = 19. The case i = refers to the weighting function 
without tomography, w| C) (^) = W r (y). One can easily observe the 
order of the polynomial hierarchy at the high distance end of the 
functions, where one after another approaches zero. 

The modified spectra C^(f) and C%(€) are shown in Fig. [4] 
and Fig. [5] respectively. The drop in amplitude is mainly an ef- 
fect of the absence of normalisation, while one can in fact observe 
slight differences in shape. However, these differences are small, 
since the polynomials only mildly depend on the multipole order 
t . Therefore, the overall shape of the spectra is still dominated by 
the zero-order spectra CS? (f) and Cf y \£), respectively. Thanks to 
the orthogonalisation these spectra now provide statistically inde- 
pendent information. In the next section we aim to combine singals 
from the galaxy distribution, the iSW-effect and the cross-spectra 
to investigate statistical bounds on cosmological parameters. 



< r w w r 0) *(O> = (fry s D ({ - {')s\ l »{0 

withS^r^octfy. 



(26) 



3.3 Properties of orthogonal polynomials 

In Fig.Q]the orthogonal polynomials are shown up to a polynomial 
order of i = 4. They show an increasing number of zero points 



4 STATISTICS 

This section aims to connect cosmic variance and statistical noise 
with the iSW-signal and its cross-correlations into a meaningful sta- 
tistical formulation. In the course of this we construct covariance 
matrices for the polynomial-weighted spectra. Statistical errors on 
cosmological parameters are estimated in the Fisher-matrix formal- 



© 2008 RAS, MNRAS 000,[T]l9] 



iSW-effect with orthogonal polynomials 5 




100 300 
multipole order t 

Figure 4. Pure galaxy-galaxy spectra S~S(f), i = 0...8, weighted with or- 
thogonal polynomials Piix), as a function of the multipole order f. Syl (■£) 
(blue) refers to the non-tomographic spectrum S yy (l). One can see a de- 
crease in amplitude for increasing multipole order I. 




100 300 
multipole order I 

Figure 5. Galaxy-iSW cross-spectra 0%{€), i = 0...8, weighted with or- 
thogonal polynomials Pi(x), as a function of the multipole order /. cf y (i) 
(blue) refers to the non-tomographic spectrum C Ty (C). 



ism. Furthermore, we investigate the signal strength of the different 
spectra and their dependence on the number of polynomials used. 



4.1 Variances of galaxy number counts 

For forecasting statistical errors, we need to derive expressions for 
the signal covariance and noise. We will start from a discrete formu- 
lation with a set of weighting coefficients w m for the counted galaxy 
number m. Clearly, the weighting coefficient w„, will depend on the 
distance of the respective galaxy. Later, we will generalise the for- 
malism to the continuous case, in which the weighting procedure 
is performed by the polynomials ptix)- The standard deviation cr ww 
of a weighted galaxy count with weighting coefficients w m is given 
by 



1 



Zm Wm £„ W„ 



/j Wm W„ S m 



(27) 



which reduces to a Poissonian result in the case of w m being either 
Oor 1: 



cr? ni , = = with n = V ' w„ . 
n *—L 



(28) 



The counted quantity n in our case is defined as the mean density of 
galaxies per steradian, for which we will substitute 7? = 40/arcmin", 
which is characteristical for the EUCLID galaxy survey. Consider- 
ing two different sets of weighting factors w m and v„, we generalise 
the standard deviation to 



1 



(29) 



which will in the continuum limit be a cross variance weighted with 
two different polynomials. For the continuum limit the transition 



dxnix)- 



(30) 



is performed which conserves the total number count n due to 
the unit normalised galaxy distribution function n(x). The discrete 
weighting sets w m and v„ are then represented by pfy) and pjix) 
so that the noise covariance in the continuous case reads 

i StcWmtoPM 

The noise term N„(i) still depends on £ since the polynomials 
are constructed for each multipole order separately. We omit the 
^-dependence of the polynomials p-,(x) for clarity. Eqn. J3 1 1 moti- 
vates the following choice of normalisation for our polynomials: 



Piix) 



Pi(x) 



J 4* nix) Pi(x) 



In this normalisation the galaxy number noise reads 



A X n (x)Pi(x)Pj(x) . 



while the galaxy spectrum can be written as 



J"XH Ay 
o X 



tlx) 



(32) 



(33) 



(34) 



The limitation in polynomial order due to increasing noise in the 
polynomials pAx) can already be illustrated: Since n(x) is a slowly 
varying function the rapid oscillations of high order polynomials 
will drive the values of the integrals j &x n 0c) Pi(X) to smaller num- 
bers and will therefore increase the noise in Piix). We point out 
that for order zero the non-tomographic case is recovered, giving 
the standard Poissonian expression for the noise iVoo = lfn and the 
integrated galaxy spectrum in the signal part Sf^iC) = S yy (£). 

While the orthogonalisation procedure leads to a diagonal 
galaxy signal covariance, the noise part will not be diagonal any 
more: t for i + j. In contrary to this method, a traditional 
binning in z would lead to a diagonal noise contribution and off- 
diagonals in the signal part. 



4.2 Fisher analysis 

In order to use both iSW signals and galaxy spectra in our Fisher 
analysis, we now define an extended data vector 



T(0 

y { Ht) 



(35) 
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The total covariance matrix, C(£) = S(£) + N(£), for these data 
vectors is block-diagonal due to the independence of the f -modes: 
Each block 



C(£) 



, v ""(£) 

j Ty\^j yy v ' 
consists of a signal part 



■"TV 

C ( ' 7) 



S(£) 



Srr(£) 



c 



where S™ (€) oc by construction, and a noise contribution 











(36) 



(37) 



(38) 



with polynomial orders < i, j ^ q. Due to uncorrelated noise in 
the CMB and the galaxy density field the noise of the cross-spectra 
vanishes. The CMB part consists of the pure iSW signal 

' %W 2 T (x)P^(k = tlx) (39) 
X 

with P^(k) = PssKxh k) 4 while the noise can be split into the pri- 
mary CMB fluctuatio ns Ccmb(Q ap d an instrumental noise term 
Cbeam(0 of PLANCK JPlanck Collaboration et al]|201 lh : 



SrAO 



T 

J() 



N„(€) = C CMB (£) + w- 1 exp (-A0 2 f) 



(40) 



with the beam size AO 

-l 



8. 77 x 1(T 4 and the squared pixel noise 

w ' = 0.2^K/rcMB (Knox 1995). The noiseless cross spectra are 
formed by one modified weighting function only: 

rxn dx 



f 

Jo 



w T (x)wf( x )P 6 t{k = eix) 



(41) 



o X 

with Ps4,(k) = Pss/ikxh) 2 - We point out that only the galaxy part 
of the signal covariance was diagonalised by our method. Conse- 
quently, the cross-spectra C®(£) are the only off-diagonal entries in 
the covariance matrix. 

The likelihood for observing these Gaussian-distributed 
modes x{£) for a given parameter set p is defined as (Tegmark, 
1997): 

£(x{£) | p) = 1 exp j 

^{2n) N detCf/) 

Defining the data matrix as Dy(€) = Xj(6)Xj(£) with (D) = C 
and using the relation, lndet(C) = trln(C) one can write the^ 2 - 
functional £, oc exp(-^ 2 /2), with help of the logarithmic likeli- 
hood L = - ln.£: 



-x T (i) C- [ (£) x'(£) \ . (42) 



X 1 = -2L = tr^ [lnC + C _I D] 

f 



(43) 



Each multipole £ provides (2( + 1) independent m-modes. If we 
interprete £, as a Bayesian probability, the local behaviour of the 
likelihood function around the point of maximum likelihood is de- 
termined by the Hesse matrix of L at this point: 

d 2 L 

(44) 



dp,! dp v 



The Fisher information matrix is then given as the expectation value 
of this quantity summed over all multipole orders £: 

F = /-*L- 

\dPpdPv 



E^f^^H" (45) 



For each £ the (2£ + l)/2 m-modes provide statistically indepen- 
dent information. In the course of our Fisher matrix calculations 
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Figure 6. Sensitivity Jtr(di, In C(l)) 2 of the Fisher matrix with respect to 
the Hubble parameter h as a function of the multipole order £, and cumu- 
lative polynomial order q. Sensitivities are shown with derivatives of all 
spectra taken into account (solid lines) in comparison to the case where 
only the cross-spectra were considered in the signal part (dashed lines). For 
the covariance the survey properties of EUCLID have been assumed. 
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Figure 7. Sensitivity ytrfflj^ In C(£)) 2 of the Fisher matrix with respect to 
the initial slope of the power spectrum n s as a function of the multipole 
order {, and cumulative polynomial order q. Sensitivities are shown with 
derivatives of all spectra taken into account (solid lines) in comparison to 
the case where only the cross-spectra were considered in the signal part 
(dashed lines). For the covariance the survey properties of EUCLID have 
been assumed. 



we will work in the limit dSjj/dp p » dNjj/dp^ and therefore ne- 
glect the noise dependence on the cosmological parameters. This 
approximation is well justified in our case. 

Next we have a look at the ratio of the sensitivities of the spec- 
tra with respect to cosmological parameters and the covariance. 
This quantity equals the contribution of a certain multipole t to 
the respective Fisher matrix diagonal element: 



8 In C(l) r 
J 



dF„, 



2^+1 At 



(46) 
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Figure 8. Sensitivity -^tr (<3n m In C(l)) 2 of the Fisher matrix with respect to 
the matter density parameter fl m as a function of the multipole order {, and 
cumulative polynomial order q. Sensitivities are shown with derivatives of 
all spectra taken into account (solid lines) in comparison to the case where 
only the cross-spectra were considered in the signal part (dashed lines). For 
the covariance the survey properties of EUCLID have been assumed. 
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Figure 9. Cumulative signal to noise ratio £ depending on the multiple- 
order £ for the survey characteristics of 2MASS (dotted), SDSS (dashed- 
dotted), EUCLID (dashed) and NVSS (solid). Shown is the improvement 
between cumulative polynomial order q = 1 (blue) and q = 8 (green). 



In Fig.[6]- Fig.[8]we show these sensitivities in solid lines for the full 
information from galaxy spectra, cross-spectra and iSW-effect in- 
cluded for the parameters h, n s and O m , respectively. At zero order 
they all exhibit a certain f -range at which the covariance is insen- 
sitive to variations of the respective cosmological parameter. Nat- 
urally, angular scales in the vicinity of this zero point do not con- 
tribute much sensitivity to the Fisher matrix. This effect is cured if 
we include all polynomials < i < q. The combination of multiple 
line of sight-weighted measurements lifts the sensitivities at these 
points continuously with increasing number of involved polynomi- 
als until the effect saturates. 

For multipole orders £ reaching higher values {£ « 3000) the 
sensitivity starts to drop rapidly. On these small scales the noise 
contribution begins to dominate and the Fisher matrix ceases to 
grow further. 

In dashed lines the sensitivities are shown if only the cross- 
spectra are included in the derivatives. Again the sensitivities grow 
with increasing cumulative polynomial order q, although in this 
case the zero order sensitivity does not suffer from any singular 
effects. Characteristic properties of the iSW-effect are recovered 
showing it to be a large scale effect due to the AT 2 proportional- 
ity originating in the Poisson equation. Above moltipoles of about 
I « 100 the information provided by the cross spectra becomes 
negligible. Clearly, the cross-spectra Fisher matrix is most sensi- 
tive to the matter density parameter C2 m , which shows the strongest 
increase in sensitivity for increasing cumulative polynomial order 
q. 

4.3 Signal to noise ratio 

A signal's power to constrain cosmological parameters is most re- 
liably quantified by the signal to noise ratio 

£ 2 = / sky V ^-^tr(c-'CClS(f>) 2 . (47) 

e 1 

Besides the physical process the signal to noise ratio strongly de- 
pends on the characteristics of the survey at hand. In the case of 



galaxy surveys the most important survey parameters are the sky 
coverage f s i y and the median redshift z mcd . In Fig [9] the signal 
to noise rati os for the survey cha racteristics of EUCL ID, 2MASS 
(z med « O.l. lAfshordi et alj|2004l) SDSS (z med » 0.5.lBielbv et al.l 
l20ld) and NVSS (z med * 1.2. iBoughn & Crittenden! 120051) are 
shown. Clearly, also the signal to noise ratio increases for higher 
polynomial orders due to the diagonal structure of the signal co- 
variance. We find an improvement of as 16% in the signal to noise 
ratio between cumulative polynomial order q = 1 and q = 8. As 
expected, at a multipole order of a few hundred the cumulative 
signal strength saturates as a result of the Poissonian k damping 
term in the iSW-effect. The actual realisation of the matter distri- 
bution in the observed universe introduces a systematic noise in 
the iSW detections known as local variance. Using the so called 
optimal method one can decrease this bias by working conditional 
on the large scale structur e data and gain 7% in signal to noise ra- 
tio IFrommert et al.l2008l) . The tomographical method presented in 
our work should be also applicable to the reconstructed large scale 
structure used in the optimal method. Therefore, a combination of 
these two methods would be sensible. 



4.4 Statistical errors 

The Cramer-Rao inequality introduces a lower limit on the 
marginalized standard deviation of the estimated cosmological pa- 
rameters. These are given by the diagonal elements of the inverse 
Fisher matrix: 

o-,>((F-\,f . (48) 

In Fig. [10] these errros are depicted for the five cosmological pa- 
rameters fl m , <Tg, h, n s and w. The plot follows the evolution of the 
errors with increasing number of included polynomials q. While for 
small polynomial orders the Cramer-Rao errors decrease rapidly 
with different characteristics for each parameter, the improvement 
slows down for higher order polynomials and assumes a character- 
istical behaviour. This behaviour can approximately be described 
by the inverse root of the polynomial order oc 1/ yTJ. A simi- 
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cumulative polynomial order q+ 1 

Figure 10. Lower limits on the marginalized statistical errors cr^ on the 
estimates of the cosmological parameters Sl m (circles), erg (squares), h 
(lozenges), « s (triangles, pointing up) and w (triangles, pointing down) de- 
rived from the Cramer-Rao inequality, as a function of the cumulative poly- 
nomial order q. The Fisher matrix was derived including the derivatives of 
all spectra C^(f). cf y (i) and C TT (t). Again, the EUCLID survey charac- 
teristics have been used. 
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cumulative polynomial order q+ 1 

Figure 11. Conditional statistical errors <x,, -C on on the estimates of the cos- 
mological parameters Cl m (circles), o"8 (squares), h (lozenges), n s (triangles, 
pointing up) and w (triangles, pointing down). The Fisher matrix was de- 
rived including the derivatives of the cross-spectra C^(f) only, EUCLID 
survey characteristics have been used. 

lar characteristic was found in the application of this meth od to the 
weak lensing shear spectra ( Scha efer & Heis enberg 201 1). Clearly, 
the cosmological parameter £l m profits most of the tomographic 
method, which was already indicated by its sensitivity improve- 
ment discussed in Section l4~2l Going to even higher orders is dif- 
ficult due to cumulative errors in the Gram-Schmidt orthogonalisa- 
tion method. 

If we are interested in how a single cosmological parameter 
can be constrained assuming that all other parameters are fixed, we 
have to study the conditional errors. These can be obtained from 
the inverse diagonal elements of the Fisher matrix 

cr„,co„ = (F^Ti . (49) 




-0.91 -0.9 -0.89 

W 



Figure 12. The 2-dimensional lcr-error ellipses for the cosmological pa- 
rameters £2 m , erg, h, n s and w from EUCLID using tomography with or- 
thogonal polynomials are shown in this plot. The lcr confidence regions 
decrease in size with increasing number of included polynomials, reach- 
ing from q = 2 (blue) to q = 8 (green). The ellipses are evaluated with a 
maximum multipole order of ^max — 3000. 

For studying the improvement provided by the cross-spectra, we 
plot the conditional errors as a function of cumulative polynomial 
order q in Fig. QT] Here, only the derivatives of the cross-spectra 
were taken into account in the Fisher matrix calculation. Again £2 m 
is subjected to the strongest improvement, its conditional error de- 
creases by as 30%. In contrast to the marginalized errors the evolu- 
tion of the conditional errors does not show a 1 / yfq behaviour but 
rather saturates at polynomial order of q « 5. 

Finally, we are interested in the 2-dimensional marginalized 
logarithmic likelihood around the fiducial model p M 

Xm \ P*-p vM )\ (F-% (F-') vr J\ Pv - PvM j ' J ' 

for which the lcr-error ellipses are depicted in Fig. [12] Starting 
with q = 2,we have combined up to nine polynomials. Besides the 
expected shrinking of the ellipses for higher numbers of included 
polynomials, it is interesting to see how the degeneracies slightly 
change their orientations in the course of tomographic improve- 
ment. This is very likely due to distance dependencies of the signal 
sensitivities. 



5 SUMMARY 

In this paper a tomographic method for measuring iSW-galaxy 
cross-spectra and galaxy spectra has been presented. It has been 
carried out by sight-weighting of the iSW-effect and the tracer den- 
sity field with specifically constructed orthogonal polynomials. 

(i) The Gram-Schmidt orthogonalisation procedure has been 
used to construct orthogonal polynomials in order to diagonalise 
the weighted galaxy signal covariance. The method projects out 
statistically independent signal contributions at the price of off- 
diagonals in the noise part. It differs from traditional tomograph- 
ical approaches, for instance from most tomographical techniques 
in weak lensing measurements, in which the noise part is diago- 
nalised. Due to cumulative numerical errors with increasing poly- 
nomial order, this method is limited to order i ss 8. 
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(ii) The improvement of the signal to noise ratios with cu- 
mulative polynomial order wa s investi gated for the galaxy surveys 



mutative poiynomial order was investigated tor tne galaxy surveys 
2MASS dAfshordi et alj|2004) . SDSS teielbv et alfeoiOh . NVSS 
(Boughn & Crittenderj|200i) and EUCLID. The signal to noise ra- 
tio for the cross-spectra only has been improved by 16% at a cumu- 
lative polynomial order of q = 8. 

(ii) A Fisher-matrix analysis was used to forcast how well cos- 
mological parameters can be constrained by different galaxy sur- 
veys, combining signals from the iSW-effect as well as from the 
tracer density field. The marginalised errors show simple inverse 
square-root behaviour with increasing number of included poly- 
nomials, which can be interpreted also as a sign of the statistical 
independent signal contributions. Conditional errors on parameters 
contrained only by the cross-spectra decrease by up to ss 30% in 
case of the matter density parameter Q m . 

(iv) While for the cross-spectra only the conditional errors 
show a saturation already at quite low number of included poly- 
nomials q « 5, it would still be worth improving this method in 
order to reach even higher orders, since marginalised errors for the 
full signal did not yet saturate at cumulative order of q = 8. 

(v) Using the wrong model in the construction of the poly- 
nomials can introduce an estimation bias on cosmological param- 
eters. This effect was thorough ly studied in a similar app roach 
for weak lensing measurements dSchaefer & Heisenberd[201lh . In 
most cases the bias was found to be small compared to the statis- 
tical errors. Since, in addition, iteration between parameter estima- 
tion and polynomial construction is able to further reduce this bias, 
a wrongly chosen cosmology appears unlikely to affect measure- 
ments. 
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